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Abstract. For finite dimensional CMV matrices the mixed inverse spectral 
problem of reconstruction the matrix by its submatrix and a part of its spec- 
trum is considered. A general rational interpolation problem which arises in 
solving the mixed inverse spectral problem is studied, and the description of 
the space of its solutions is given. We apply the developed technique to give 
sufficient conditions for the uniqueness of the solution of the mixed inverse 
spectral problem. 



1. Introduction 

The theory of CMV matrices, rapidly developing in the recent years ([UlSllSjj see 
also the expositions in [131 [3] and the references therein), has a strong background 
in the theory of Jacobi matrices. The similarity between the Jacobi and CMV 
matrices not only provides a general concern of investigation, but also permits 
sometimes to predict the most probable answers. Some of the facts known for the 
Jacobi matrices can be easily extended to the CMV case, some require considerable 
adjustment. In turn, the spectral theory of the Jacobi matrices is paralleled by the 
spectral theory for the Sturm-Liouville differential operators, which forms another 
important front of study in this area of research. 

Mixed inverse spectral problems (MISP) are of special interest in the inverse 
spectral theory. For this kind of problems one reconstructs a differential or differ- 
ence operator by a part of its potential and some additional spectral data. Started 
by Hochstadt and Lieberman [7] for the Sturm-Liouville operators, these problems 
were extended and refined for many other cases (see [4] for the references). Com- 
pared to the "ordinary" inverse spectral problems where the whole potential is to 
be reconstructed, for MISP one needs to know "less" spectral data. 

In what follows, we refer to several key studies of the development of MISP. Let 
J be an X iV three-diagonal matrix of the form 



/ai 6i 
bi a2 62 
62 as 63 



bN-i aN J 



a„, 5„ e M, bn > 0. 
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Consider the a's and the b's as a smgle sequence ai, 5i, 02, 62, • ■ •, that is, 

C2n-i = a«, C2„ = 6„, n e N. 

In [S] Hochstadt proved the discrete version of the Hochstadt-Lieberman theorem: 

Theorem. Let n E N. Suppose that cn+i, ■ • ■ , C2n~i o,fe known, as well as the 
eigenvalues Ai, . . . , Aat of J . Then ci, . . . ,cn are uniquely determined. 

It is important to remark that in certain compHcated physical systems it is 
not always possible to know the entire spectrum. A natural question arises if it 
is possible to reconstruct the Jacobi matrix when we know more than a half of 
the potential, but less than the whole spectrum. The positive answer is given by 
Gesztesy and Simon in [F. 

Theorem. Suppose that I < j < N and Cj+i, . . . ,C2n-i are known, as well as 
(any) j of the eigenvalues. Then ci, . . . , Cj are uniquely determined. 

The goal of this paper is to study the MISP for CMV matrices. Although the 
algorithm of solving the problem is similar to that of the Jacobi case, essential 
difference arises when the uniqueness of the reconstruction is concerned. In the 
Jacobi case the MISP is reduced to the interpolation of a rational function (specif- 
ically, the Weyl function of the unknown submatrix) by its values in the known 
eigenvalues of the whole matrix. The degree of the rational function (i.e., the sum 
of the degrees of the numerator and denominator) corresponds to the number of 
the interpolation points, and both numerator and denominator are monic. So, the 
uniqueness of the interpolating rational function drops out immediately from the 
simple fact that a polynomial of degree k with fc + 1 zeros is identically zero, and 
thus we prove the uniqueness of the reconstructed Jacobi matrix. However, in the 
CMV case the degree of the interpolating rational function (the Weyl function of 
the "reduced CMV matrix") is greater by 1 than the number of the interpolation 
points, and the lacking "piece of information" is given by a restricting condition on 
the free term of the numerator to be 1 . Here the trivial consideration of the Jacobi 
case fails, and a theory is required of how to find the interpolating rational function 
with such restriction. 

So, the work consists of two parts. In Section 2 we give an approach to the 
rational interpolation theory adapted for solution of the interpolation problem re- 
lated to the MISP. The MISP itself is studied in Section 3. Note that, although in 
the Jacobi case we always have uniqueness in the MISP (provided the number of 
lacking entries agrees with the number of the given parameters), in the CMV case 
a degenerated case is possible where the MISP has infinitely many solutions. The 
main results of this paper are Theorems 12.121 with a description of the solutions of 
the rational interpolation problem, and 13. 31 with a sufficient condition for the MISP 
to have a unique solution. 

2. Two-dimensional vector-polynomials and rational interpolation 

We start with the following interpolation problem: given points zi , Z2 , . . . , G C 
and numbers wi , 012 , • ■ • , w„ G C find a "nice" description of all rational functions 
p(i)/p(2)^ where P^^^ and P*^^^ are polynomials with complex coefficients, for which 

(2-1) 75(5)7r(='^^-' = 
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{ujj — oo means that the rational function must have a pole at the point zj). By 
a "nice" description we understand a description whose form allow us to "control" 
the degrees of P^^\z), j = 1,2, and some more parameters such as their leading 
coefficients or the free term, required in the concrete application od the rational 
interpolation problem (we will return later on to this topic). 

It is advisable to reformulate this problem as a linear problem in the space of 
two-dimensional polynomial vector- functions (vector-polynomials) . Define numbers 
a^j''' and a^^'' by the following rule 

\ aj^i := 0, aj,2 := 1, if i^j — oo. 

Then ((2?T|) implies 

(2.3) a,-iP«(zj) + «j,2P^'H^,) = 0, j = l,...,n, 

where -I- \ctj^2\ > 0. 

Conversely, if \P^^Hzj)\ + iP^-^^Zj)] > 0, j = 1, . . . , n, then (EJ]) implies (HH) 
with 

(2.4) ujj-.^-^eC, j^l,...,n. 

However, if \P'^^\zj) \ + \P'^^\zj)\ = for some j, then (|2.ip does not make sense 
at the point Zj. 

Thus, we reformulated the initial interpolation problem in the following way: 
find all the pairs of polynomials P^^-* and P'-^\ for which (|2.3p holds. The second 
problem is "almost equivalent" to the first one: if the rational function pW J pi"^) 
is a solution of (|2.ip . then the pair of polynomials P'-^-' and P'^^-' is a solution of 
p.3p with aj^i and aj,2 defined in (|2.2p . If P^^^ and P'^^^ are polynomials such that 
|P(i)(zj)| \P^'^\zj)\ > solving ([23]), then the rational function pW/pC^) is a 
solution of ()2.ip with w^- defined in (12. 4p . 

In what follows we refer to interpolation problem (|2.3p as the problem (/„), and 
denote by Poo(-^n) the set of all its solutions. Some of the results desribed below for 
the polynomial vector- functions are taken from [lOj . [llj where the interpolation 
problems also appears as well as its main objects like generators, etc. 



2.1. The space of vector-polynomials. To describe the class of solutions of 
we introduce the space 

/ P(i)(z) \ 1 
p{z) = I p(2)| j ) ' ' * = li2, are complex polynomials | . 



This is a linear space with the standard operations. The zero element on this space 





is = ( ^ ) . We point out that Poo is also a module over the ring of polynomials: 



Sp 



^P(I) 

5P(2) 



for any polynomial S. 
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p(l) 
p(2) 



^ is the 



Definition 2.1. The height of the vector-polynomial p = 
number 

r 2 deg P(i) , deg P(i) > deg p(2) ^ 
^^•^^ ^^^^ - \ 2dcgP(2) + 1^ degP(i) < degP(2), 

As usual, degO = — cxi, so we put h{0) := —oo. 

It is obvious that 

(2.6) h{p) = max{2 deg P^^^ , 2 deg P^^^ + 1} 
and 

(2.7) /i(S'p) = /i(p) + 2degS'. 

The degrees of the components of the vector-polynomials can be written down in 
the following table: 



height p 





1 


2 


3 


4 




2k 


2k + 1 




degP(i) 


= 


< 


= 1 


< 1 


= 2 




= k 


< k 




degP(2) 


— oo 


= 


< 


= 1 


< 1 




<k-l 


= k 





The following proposition demonstrates that the notion of the height is a natural 
extension of the degree of polynomials. 

Proposition 2.2. If h{p) ^ h{q), then Va,5 e C 

h{ap + bq) = iiiax{h{p), h(q)) 

U h{p) = ^(9) =^ then 

(1) y a,beC h{ap + bq) < n 

(2) 3c e C: h{p + cq) <n-l. 

Proof. We prove (2), the rest is plain. If n — 2k, then k = degP'"'^' > degP(^\ 
k = degQ(i) > degQ(2). So, there exists c G C such that deg(p(i) +cQ(i)) < fc-1. 
Also, since degP^^^ < fc- 1 and degQ^^^ < A:- 1, we see that deg(p(2) +cQ^'^'>) < 
k-1. Then, by ([^ 

h{p + eg) < 2fc - 1 = 71 - 1. 
If n = 2A: + 1, then k = degP^^) > degP(i), k = degQ^^) > degQ(i). So, 3c e C: 
deg(p(2) +cQ(2)) < fc - 1, deg(p(i) +cQ(i)) < fc. Then, by (EH) 

h{p + cq) <2k^n-l. 



□ 



Consider the following basic system of vectors in 



(2.8) e2k{z) := j , 62^+1 (z) := 

It is clear that h{en) — n for all n. 
Proposition 2.3. {e„}„>o is a basis in II 

pGPoo, h(jp)=m = 

and this representation is unique. 



k - 



3; i.e., 



/(:eZ+ = {0,1,...} 



fe=0 
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Proof. We use induction on m: 

1) The cases m — 0,1 are checked immediately 

2) Let m = 2fc + 1, then k = degP^^) > degP(^\ i.e., P^^^z) = az'' + Q^^\z), 
a^O, degQ(2)(2) < - 1. Then 

q{z) := p{z) - ae2k+iiz) = (^q(2)|^| 

where Q'^^^z) = P^^^z) and degQ^^) < k. 

So, h{q) < 2k — m ^ 1, and we can apply the inductive hypothesis. The repre- 
sentation q{z) = J2T=o Cfcefc(-z) is unique, so the representation 



p{z) = aem{z) + ^ Cfeefe(2 



k=0 



is also unique. 

3) Let m = 2k. ThenP(i)(z) = bz'' + Q^^'> (z) , b^O anddegQ^^^ < k-1. Define 
q{z) := p{z) - be2k = (^gja) • Here deg Q^^) < fc - 1, so h{q) < 2fc - 1 = to - 1 
and we can again apply the inductive hypothesis. □ 

The latter proposition can be extended in a natural way. 

Proposition 2.4. Let {gniz)}n>o be an arbitrary sequence of vector-polynomials, 
such that 



h{gn) = n, 



n E 



Then {(7„}„>o is a basis in f 
Proof. By Proposition 12.31 



k=0 

or in a vector-matrix form 

/ eo 







\9m) 





Clo cii 



\CmO 



where the matrix of {cm.k) is triangular. So, 







Vem^ 





ClO Cii 



50 



^0. 



Since {e„}„>o is a basis, then so is {5ri}n>o- 



□ 
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2.2. Transforms in Poo- A 2 x 2 matrix A = [ , ) defines a transform of tlie 



vector-polynomials : 



c d 



A f \ - A fP^^Kz)\ _ faP^^Hz)+bP^^Hz)\ 

Proposition 2.5. (1) For arbitrary matrix A and p G Poo 

h{Ap) < h{p) + 1. 

(2) // A is upper-triangular, then for arbitrary p £ Pqo 

h{Ap) < h{p). 

(3) // A is lower-triangular, then 

Hp) <2k-i-l=^ h{Ap) < 2fc 4- 1. 

Proof. (1) By (HU 

h{Ap) = max (2 dcg(aP(^) + foP^^)), 2 deg(cP(i) + dP^^)) + l) , 



2deg(aP(^) +fep(2)) < 2max(degP(^\ degP'^^) = max(2degP(^\ 2degP(2)) < h{p), 
2deg(cP(i) +dP'-^^) < max(2degP(i\ 2degP(2)) < h{p), 

2deg(cP(^) -I- dp(2^) -f 1 < h{p) + 1. 

(2) If c = 0, then 2deg(dP(2)) + 1 < 2degP(2) + 1 < h{P). 

(3) Let now h{p) < 2k + 1 and A = ("^^ |j ) . Then 



c 



. , ,( aP(^\z) \ 

^P^^>- l^cPW(z)+dP(2)(z)j 



and we have by assumption 2degP(i) < 2fc + 1, 2degP(2) < 2k. So, 

deg(cP(^) + dP'-^^) <k=^ h{Ap) <2k + l, 
as claimed. □ 

Later on we will use the following property of the height, which is a simple 

, , /P(i'\ 

consequence of (|2.6p : for p = I ^ ^00 we have 



(2.9) 



Vp(2) ; 

h{p) <2k+l^h (^^' 'pi^^^^'^^ <2k + 2, 
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2.3. The generators of interpolation problem. It is clear that solutions of 
(|2.3|) form a module over the ring of polynomials in Poo, i.e., if r and q are solutions 
of (12. 3p . then so is Sr + Tq for arbitrary polynomials S and T. The goal of this 
subsection is to show that this module has exactly two generators, and to study 
their properties. 

Recall that Poo(-^n) is the set of all solutions of (|2.3p . Set 

/i(/„) :=min{/i(g) : gePooUn), g^O}, 
which we call the height of the interpolation problem. 

Definition 2.6. We say that r —Cz Poo(^n) is a minimal generator of 1^2. 3\i . if 

h{r) = h{I,,). 

Proposition 2.7. The minimal generator of {!„) is unique up to a constant factor. 



Proof. Let ri and r2 be two minimal generators. By Proposition 12. 2| 3a £ C: 
h{ari + < h{In) — 1. But ari + br2 is a solution of (|2.3p . Since ri and r2 are 
minimal non-trivial solutions of (|2.3p . we conclude ari + br2 =0. □ 



A trivial (nonzero) solution of ((O)) P'^) = P^^) ^ Hjl^-^i) 

provides the bound 
h{-^n) < 2ri + 1. It turns out that this bound can be improved immensely. 



Theorem 2.8. /i(/„) < 



n. 



Proof. The following non-negative matrices of rank 1 play a key role in our 
consideration: 

o-fc = - I 12 = - ] {ak,i, ak,2), fc=l,...,n. 

\ak,2ak,i \ak,2\ ) \(^k,2 J 

It is clear that 

\ak,iP^^\zk) + ak^2P^^\zk)\^ = p*{zk)crkp{zk), 
p*iz) := (PWIT), W^iD), p{z) (^^I^Jl^j 



so the problem ()2.3|) is equivalent to 

(2.10) p*{zk)(Tkp{zk)=Q, k = l,...,n. 

We proceed by induction on n. 

1. For n = 1 we have a non-trivial solution 

2. Suppose that we have already proved the result for n, and we want to prove 
it for n + 1. The forthcoming construction depends on whether n is odd or even. 

Let n — 2k + 1, and consider the problem {In+i) (|2.10[) with n + 1 data. If 

aj^i = for all J = 1, . . . , n + 1, then the vector-polynomial p = , h{p) — 0, 

is a solution of (/«), and we are done. So, suppose without loss of generality that 
an+1,1 9^ (otherwise enumerate the points zi, . . . , Zn+i). The upper-triangular 
matrix 

^(Oin+l.l)^^ 2(a„+l,l)^-^^ 

1 



no := 



8 
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satisfies 

'1 



0/ ' 



(2.11) rigCTn+iilo 

and 

(2.12) n*a,no^(l^'''f ^^^i^; J = l,---,n, 



■^7 ■ ■ ■ 5 



with some numbers (3j,i, f3j,2, J — 1, ■ • • , fi- Put 

7i,i (zn+i - Zj)PjX; li,^ ■= Pja, j = 1, 

and consider an auxiliary interpolation problem (/„): 

Jj,iP^'\z,) + l,.2P^^\z,) = 0, J = 1, . . . , n. 

By the induction hypothesis, there exists a solution q e Poo, h{q) < n = 2fc + 1, of 
this problem: 

g*(zj)CTj g(zj) = 0, j = l,...,n, 

where 



Define a vector-polynomial 



(2.13) '^W-U.)y)-"o, 



For r we have 

r*{zj)(Tjr{zj) = q*{zj)aj q{zj) = 0, j = 1, . . . ,n, 
and for j ~ n + 1 

r*(z„+i)CT„+i r(z„+i) = g*(z„+i) "j*^ ^jj "j*^ (?(2;«+i) = 0, 

by (|2.11[) . So, r is a solution of (/„+i). 

Since h(q) < 2fc + 1, then by the upper inequality in (|2.9p we have 



< 2A: + 2 = n + 1, 



'{Zn+l-z)Q(^HzT 

so, by Proposition 12. 5[ h{r) < n + 1, as needed. 

Let n = 2k, and (/n+i) (|2.10p be the interpolation problem with n + 1 data. If 

Q!j,2 ~ for all J = 1, . . . , 71 + 1, then the vector-polynomial p — ^'j'^ , h{p) = 1, 

is a solution of (/«), and we are done. So, suppose as above, that Q!„+i.2 7^ 0. The 
lower-triangular matrix 

^ ' V~""+14('^"+1.2) ("„+l,2) V ' 

satisfies 

(2.14) f^ta„+ini = 
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and 

(2.15) = (£C M^)^ ^ = '■■■■■"■ 

Put 

^j,2 '■= {zn+1 ^ Zj)5j,2', Aj^i :— 5j^i, j = 1, . . . ,n, 
and consider an auxiliary interpolation problem (/«): 

A,- iP(i)(zj) + Xj,2P^-^\zj) =0, j = 1, . . . ,n, 

By the induction hypothesis, there exists a solution q G Poo, h{q) < n = 2k, of this 
problem: 

q*{zj)ajq{zj) ^0, j = 

where 

Define a vector-polynomial 
For r we have 

r*{zj) Gj r{zj) = q*{zj) aj q{zj) = 0, j = 1, . . . , n, 
and for j — n + 1 

r*(z„+i) cr„+i r(z„+i) = g*(2;„+i) Q = 0, 

by ()2.14|) . So, r is a solution of 

Since h{q) < 2k, then by the lower inequality in ()2.9p we have 



QW(z) 
(z„+i-z)Q(2)(^) 



< 2A: + 1 = n + 1, 



so, by (3), Proposition 12.51 {ili is lower-triangular), h{r) < 2k + 1 — n + 1. The 
proof is complete. □ 



If r is a minimal generator of (/«), then Sr G Poo(^n) for any polynomial 5*. So, 
the question arises naturally whether Poo(^n) — {Sr}, S a polynomial. The answer 
is negative: it turns out that the module of solutions of (/„) has exactly one more 
generator. Denote F'^{In) = Foo{In)\{Sr}. It is shown in Theorem 12.101 below 
that this set is nonempty, so the following definition makes sense. 

Definition 2.9. We say that q G Poo(^ri) *s a second generator of (In), if 

h{q)^min{h{p), peP'^(/„)}. 

Theorem 2.10. The set P^(/„) is nonempty. Furthermore, the height of any 
second generator is h{q) = 2ri -I- 1 — h{In). 



10 
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Proof. We show that there is a solution of (/„) of the height < 2n + 1 — h{In), 
which is not of the form Sr. 

Let h(In) — k < n. Pick arbitrary different numbers Zn+i, • ■ ■ , ^2n+i-fe G C 
distinct from zi, . . . , z„. Take numbers a^.i, 01^.27 j = n + 1, n + 2, . . . , 2n + 1 — /c 
in such a way that 

(2.17) aj-ii?(i)(zj) + aj-2i?^^H^j) 7^ 0, j = n + 1, n + 2, . . . , 2n + 1 - fc, 

/ \ . . . 

where ^ = ( ^^(2) j minimal generator of (/«)• Consider the interpolation 

problem 

(2.18) aj-iP(i)(zj) + aj-2P^^Hzi) =0, j = 1, 2, . . . , 2n + 1 - fc. 

By Theorem 12.81 there exists a nonzero vector-polynomial p = I ^^2) j 7^ 0> 

which solves this problem, and h{p) < 2n + 1 — k. 

Suppose that p = Sr for a polynomial S. Then, by (|2.18p . 

S{zj) [ajui?(i)(zj) + aj,2R''^\zj)] =0, j = n + 1, n + 2, . . . , 2n + 1 - /c, 

so by ((2T7l) . 

SiZn+l) = S'(z«+2) = . . . = S'(z2n+l-fe) 0. 

Since S'(z) ^ 0, we see that deg S{z) >n+l-k, and by ([2^ 

/i(p) = /i(S'r) = /i(r) + 2degS' > h{r) + 2{n+l -k) = 2n + 2~k, 

which leads to contradiction with h{p) < 2n + 1 — k. So p G ¥'^{In) and h{p) < 
2n + 1 — fc, as claimed. 

Let now q be any second generator, so h(q) < 2n + I — k. We prove next 
h{q) >2n+l-k. 
We have 

a,,ii?W(z,)+«,,2i?(2H^j-) = 0, 
a,,iQ(^Hz,)+a,,2Q^^Hz,)^0, ^ ^,---,n, 

so 

(2.19) det (^Hl-^^l q'^\z!)) = ^^'^(-.OQ^'H-,) - R^'Hz,)Q('\z,) = 0, 
j — 1, . . . ,n. Suppose that h{q) < 2n — k, which implies by (|2.6p 

degQ 

If h{r) = k is even, then 



k k 
degQ(i)<n--, degg(2) < n- - - 1. 



degi?« = ^, degi?(2)<^-2, 

and so 

deg(i?«g(2)_^(2)Q(i)) 

If fc is odd, then 

degi?(2) = ^, degi?W<^, 

and again 

deg(i?(i)g(2)-i?(^)Q(i)) <n-l. 
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It follows now from (I2.19P that 

(2.20) (z)Q(') (z) - (^)Q(i) (^) = 0. 

Let T be the greatest common divisor of R'^^^ and so (|2.20p turns into 

X(i)TQ(2) = x(2)tQ(i) , with relatively prime X^^) /T and X^^) i?(2) /T. 

Hence S'^^^ Qt^V^^^^ is a polynomial, and so is S'^^) Q(2)/^(2) js^^^ 
X(i)TQ(2) ^ x(2)TQ(i) implies ^(i) = ^(2) ^ ^ and 

(2.21) qW^:^.5, z=l,2. 

Note that all roots of T are among the nodes of interpolation. Indeed, if T{w) ~ 
and w ^ {zi, . . . , z„}, then 



r{z) = 



2)j ePoo(/n) 



and /i(r) < /i(r), which is impossible since r is a minimal generator of {In)- Hence, 

W ^ Zl. 

It remains only to show that S{zi) — 0. Assume that S{zi) ^ 0. Then 



and (P?^ imply 



So 



, i?W(z) R^^\z). ^ 
1™ {"'.I 7 + "',2 7 ^ } = 0, 

z-*zi l^z — Zl)"' [Z — Zl)"-' 

where T = {z - zi)"'f, f(zi) 0. Since = {z - zi)"'3'\ « = 1, 2, we see that 

r{z) = -. ^-r(^) ePoo(/n) 

[Z — Zl)"' 

and clearly h{r) < /i(r), which again leads to contradiction with r being the minimal 
generator. 

Finally, since all the roots of T are among {zi, . . . , (k} and S{zi) = . . . — S{zk) = 
0, P = S/T is a polynomial, and by ^^TU^ Q^"^ = PR('\ which contradicts to 
q G P'oq(/„). The proof is complete. □ 

Remark. In fact we have proven that each solution q G P'^(/„) with h{q) < 
2n + 1 — h{In) is a second generator. It is also not hard to see that each solution 
q G P(/„) with h{q) = 2n + 1 — /i(/„) is a second generator. 

Theorem 2.11. Each solution of the problem (/„) has the form 

(2.22) p{z)^S{z)r(z)+T{z)q{z), 

where r and q are the minimal and second generators of (In), respectively, and S 
and T are polynomials. Conversely, each vector-polynomial of the form l\2.22^ with 
arbitrary polynomials S and T belongs to Poo(^n)- 
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Proof. We only prove the first statement. Let — k. Consider a system of 

vector-polynomials {/j"}j">o defined as follows: 

fk — hn+l-k ~ Q, 

fj = ^j^ j = 0, 1, . . . ,/c - 1; fk+2j-i = ek+2j-i, j = 1,2, . . . ,n - k, 

fk+2j = z^r, f2n+i-k+2j = z^q; j e N, 



Bi are in ()2.8|) . It is easy to check that h{fj) — j for all j, so by Proposition 12.41 
this system is a basis in Poo, and in particular each p S Poo(^n) admits a unique 
representation in the form 



fc— 1 n—l—k 



p(z) = ^ OiBi + ^ biek+i+2i + Si{z)r{z) + Ti{z)q{z). 



i=Q 

fe-1 )l-l-fc 



Since p,r,q £ Poo(^ri) then Oje,; + biek+i+2i is also solution of (/„). But 

i=a 

its height is less then 2n+ 1 — fc, so, according to the definition of a second generator, 

fe — 1 n — 1 — fe 

^0^6^+ ^ biek+i+2i ^ S{z)r{z). 

1=0 

Hence p is of the form ({2:22)1 with S = Si+S,T ^Ti. □ 

Thus, the following theorem holds for interpolation problem (I2.ip : 
Theorem 2.12. Each solution of problem (Ig.ill has the form 

S{z)R<^^\z)+T{z)Q<^^\z) 



S{z)Ri-^){z) + T{z)Q(^){z) 



(2.23) 



where r — ^^(2)J '^"■'^ 9 ~ yQ^'^-'/ '^''^ minimal and second generators of {In), 
and S and T are polynomials. Conversely, if r and q are the minimal and second 
generators of (In), o,nd S and T are such polynomials that the numerator and the 
denominator in ^2. 23^ have no common roots, then i2.23^ is a solution of ^2.1]i . 



Remark. Roughly speaking, the task of giving a description for the set of the 
solutions of rational interpolation problem ()2.ip . is obvious. Let 'j^j2j ^'^ '^'"^H ratio- 

nal function solving problem (|2.ip . Then all the functions of the type 



i?(2) Q(2)' 
... 

where q(2J arbitrary rational function vanishing in the nodes of interpola- 

tion, will be all the solutions of the rational interpolation problem. However, such 
a description does not permit us to predict the degrees of the numerator and de- 
nominator as well as other properties needed in applications. So, we cannot obtain 
in this way a rational function which solves the interpolation problem and has the 
prescribed properties. For example, in the interpolation problem appearing in the 
next section we will need a rational function with monic numerator and denomina- 
tor of concrete degrees, such that the free term of the numerator equals 1. Thus, 
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more elaborated results are required. Certainly, we do not think that the two de- 
scriptions for the solutions of the interpolation problem, mentioned above, are the 
only possible. 

3. Reduction of MISP to Rational Interpolation 

For the definitions, notations and basic properties of finite CMV matrices see, 
for example, [HI [6]. We will add some more to the list. 

Let C = C(ao, . . . , ctn-2] P) be a finite CMV matrix with Verblunsky's parame- 
ters (ao, . . . , Q!„^2; /3) and the system of the Szego polynomials {$o, ■ • • , ^n-i', ^n}- 
They satisfy the famous Szego recurrence relations 



(3.1) 



$fc(z) = z$fe-i(z) - afc_i$fe_i(2:). 



l,2,...,n-l, $0 = 1, 



<f„(z) = z$„_i(z)-/3$;-iW- 
As is known, 

71 

a spectrum of C, (j Q, j i. Put 

m— 1 



(3.2) 



3=0 



2N-1/2 



0,1,. 



and for appropriate values of the indices define 

(3.3) X2kiz) -.^ Z^'^K2k^2k{z), X2k+l{z) ■■ 



-k-1 



Kq = 1, 



'«2fc+l$2fe+l(^) 



where by definition (filiz) = z''(pk{l/z). 

For each eigenvalue Q the following equality 

(3.4) CX, = QX„ X, = [x^{Q, • ■ • , Xn-i{Q]\ 

gives (along with (|3.3p ) an explicit expression for the eigenvectors of C in terms of 
the Szego polynomials and Verblunsky parameters. p.4p is proved in |121 Lemma 
4.3.14], for infinite CMV matrices. For finite matrices the argument is similar. As 
a matter of fact, the following more precise result holds. 

Proposition 3.1. For z e C\{0} and X{z) = [xq{z), . . . ,Xn~i{z)]* the equality 
holds 

(z-C)X = z-["/21k„_i$„(z)«„, t;„eC", ||«„|| = 1. 

Due to the sieving procedure we will assume without loss of generality that n 
is an even number: n = 21. Throughout the rest of the paper we assume that 
the last Verblunsky coefficient is known, and for simplicity put /3 = 1. The CM 
factorization now takes the form C(ao, . . . , ct2i-2, 1) — with 
(3.5) 

^^(ao) \ e{ai) 

M = 



C 



(0(21-2)/ 



{(^21-3) 



1/ 



e(a,) 



> 0. 
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Put 



(3.6) 



U = 



7 



J = 



1 

1 



the orthogonal 21 x 21 matrix, and consider the reflection of C 
(3.7) Cr UCU = UCU ■ UMU = Cr ■ Mr- 

It is clear from (|3?5)) - ((3J)) . that 

-a2/-2) 



(3.8) 



Cr = 



Mr 



0{-ao)J 



/I 



V 



'0!2l-3) 



0{-ai) 



(3.9) Cr = C(Ao, . . . , A2/-2; 1), Afe := -a2/-2-fe, k = 0,1, ... ,21 ~ 2, 

is also a CMV matrix corresponding to the "reversed" Verblunsky parameters. We 
denote the Szego polynomials for Cr by {Aq, . . . , A„„i; A„}. 
Obviously, Y.(Cr) = S(C) = {Q}" and 

CrY^ = QYj, = [y„{Q, Vn^QY, 

where yk are in (|3.3p for the matrix C^. On the other hand, by (|3.4p and (|3.7p 

= = [xn-i{Q), ■ • ■ , a:o(0)]*- 

Since the spectrum is simple, the vectors Yj and Xj are proportional: 



CjXj-, yfc_i(Cj) = Cja;„_fc((j); 



1,2,. 



or 

(3.10) 



2/fc-i(Ci) _ x„-kiCj) 



VkiCj) Xn-k-liCj)' 



fc = 1,2, . . . ,71 - 1. 



Under the mixed inverse spectral problem (MISP) we mean the reconstrucion 
of a CMV matrix C — C{ao, . . . , ctn-2', 1), or equivalently, of a set of Verblunsky 
parameters ao, . . . , a„_2 G IC, when a part {Cj}JLi of its spectrum and a part of 
the system (ao, . . . , an-2) are known. 

Here is the simplest problem of this type. Assume that we know (ao, ■ • ■ , ctn-s) 
as well as two eigenvalues Ci C'l, and an-2 is to be found so that ^1^2 G ^(C). 
Once $„-2 is known, we apply the Szego recurrences to obtain 



so 

(3.11) 



$„(z) = Z{z + a„_2)*n-2(^) - {zan^2 + l)*,*_2. 

A + a„_2 



biQ=T„ 6(A) 



1 + Xan 



J = 1,2, 
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The question is whether a„_2 is uniquely determined from the interpolation prob- 
lem p. lip . An elementary analysis of ()3.1ip shows that it has a unique solution as 
long as nCi ^ T2C2, that is, 

K-2{Cl) ^;-2(C2) 
^ ^ ^ <i>n-2(Cl) $«-2(C2) ' 

it has infinitely many solutions if T2 = — Ci ^^'^ '''i — ~C2, or 

^;-2(Ci) ^ ^;,2(C2) ^ _ 

$n-2(Cl) <i>n-2(C2) ' 

and it has no solutions at all, if = T2C2, but ti 7^ — C2, ''2 7^ ~Ci- It is not hard 
to check that each situation may occur for interpolation problem p. lip . However, 
if the existence of CMV matrix C with Ci_2 S S(C) is supposed, the existence of the 
solution of problem (13. lip is guaranteed and the problem of finding q;„_2 may have 
either unique or infinitely many solutions (see example 1 below). 

Since the Blaschke product $„_2/*&^-2 of order n — 2 cannot take the same value 
on the n-point set S(C), there always exists such a pair ^ ^2 in S(C), that p.l2p 
holds, so a„_2 is uniquely determined. 

The general MISP for CMV matrices we study here looks as follows. Let n = 21 
be even. Given first n—m—1 Verblunsky parameters oq, . . . , an-m-2, and 2m eigen- 
values Ci, . . . , C2m, 1 < m < n/2 = I, find the rest m parameters q;„_„i_i, . . . , a„_2 
and thereby restore the whole matrix C. Our main result provides the conditions 
for this problem to have a unique solution. 

Consider a pair of CMV matrices with the ''known" parameters C(q;o, ■ ■ ■ , ctn— m— 3i 
and C(ao, . . . , ctn-m-2] 1) and the systems of the Szego polynomials 

{$0,...,$ n — rn — I ;$n-^m~l}, {$0,...,$ n — m — 1 

respectively. By the Szego recurrences (13. ip 

'^n-m-l{z) = Z^n-m-2{z) ~ $:_,„_2(^), 
*ri-m-l(z) = Z^n-7n-2{z) ~ a„_„, -2 *Jl-m-2 (-2) i 

SO 

(3.13) $„_,„_i(z) - $„_,„_i(z) (1 - an-„i-2Wn-m-2- 

Similarly, for the pair C(Ao, . . . , \m-i \ 1) and C(Ao, . . . , Xm \ 1) of "unknown" CMV 
matrices, \j from (|3.9p with the Szego polynomials {Aq, . . . , A™; Am+i} and 
{Ao, . . . , A,„+i; A„i+2}, respectively, one has 

(3.14) A™+i(z) - A„+i(z) = (1 - A„0A;,. 
Now write p.lOp with fc = m + 1: 

(3.15) -^^^ = ^"—^j^^j , j = l,...,n, 

ym+l(y) a:^n-m-2(0) 



^ 2m "real" parameters are given to find m "complex" ones. 
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and observe that the right hand side of (|3.15|) is known for j ~ 1,2, . . . , 2m. Indeed, 
let, e.g., m be odd (for even m the calculation is the same). Then by p.3p 

Xn-m-2{z) = Z 5 Kn_,„_2$,* _m-2 (^) 5 Xn-m-l{z) = Z 2 Kn-m-l^n-rn-l{z) , 

so in view of p.l3p 

•^n— m — l(Ci) ^n— rn— 1 ^n— m— l(C:/) 



Xn-m-2{Cj) Kn-m-2 *,*-m-2(0) 
= (1 - \an-m-2\'^) 



2^-1/2 + (1 - «»-m-2)^r*-m-2(0) 



Pn-m-2 



^n-m-l{Cj) -. I ^ /'I l„ |2\l/2 

— - + 1 - a„-™-2, > , = (1 - lad ) • 

^n-m-2\^j ) I 



In the same way 



and with 



we have 



n(i-iA,n"^/^= n (1 

J— n — 771 — 1 



|2n-1/2 



Pm,r 



ym+i(0) 



Am+l(Cj) -, I _n IX |2^-l/2 _ „ 

+ i + an-m-2 } , Pm,r — U ^ |^m| J — Pn-m-2 



Using p.lSp . we end up with the following equalities for j = 1,2,..., 2m 



-m-2| 



|2 



(3.16) = -1 - an-m-2 + ^ 

As the last step, we express the ratios in terms of the Weyl functions (more 
precisely, their reciprocals) 

N 1 $„-™-i(z) I ^ 1 Km+ii.z) 

W{z) = — - = , Wr[z) = = — 

W{z) $„_™_2(z) Wr{z) Km{z) 

of the known C{aQ, . . . , an^m~z] 1) and unknown C(A, . . . , Am-i; 1), respectively. 
Indeed, for \z\ — 1 

n — 7n — 1 n—rn 

$„_„_i(z)= II {z~z,), = $„_™_i(0) = (-l)"-™-i [| = 

71 — 7n — l 



so* 1 fiT'l = -<r~"+"'+^<I> ifiTO 



$„_„_l(C,-) _ g- <^n-m-l{Q 
K-m-2iQ C""""'<i>n-™-2(0) 
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and, similarly, 

Am+l(0 



Finally, we come to the following interpolation problem for the Weyl function of 
the "unknown" CMV matrix C(Ao, . . . , Xm-i', 1) 
(3.17) 

„, f.,. (1 + an-m-2) {1 - " CjW {Q)} " (1 " \an-^-2\^) _ 

^r((,jj - (,j J- , y. UJj , 

|1 - an-m-2 - 

j ^1,2,.. .,2m, or 

(3.18) P«(0)-C.,P(2)(0) = 0, j = l,...,2m, 

with LUj defined in (|3.17p . which we have denoted by (/2m) in the previous section. 
Now LUj ^ 00 since all zeros of A,„ are in the open unit disk D. The above argument 
shows that (|3.18[) has a nontrivial solution 



A(2) ; " V A 

and 

(3.19) /i(A) = 2m + 2 

Proposition 3.2. For the problem p.lSp /i(/2m) > 2m — 1. 

Proof. Let r be the minimal generator of (|3.18p . and suppose that h(r) < 2m — 2. 
By Theorem 12 . 101 for the second generator q one has h{q) > 2m + 3. It follows now 
from Theorem and ([XTO)) that X = Sr with degS" > 2, so 

A„+i(z) = S{z)R<^'\z), A™(z) = S{z)R^^Hz), 

which is impossible, for Am+i and Am have no common zeros. □ 
There is some more information available about the solution A. Specifically, 

(3.20) deg A„i+i = m + 1, deg A„ = m 
and 

(3.21) A„+i(0) = -1. 



In view of Theorem 12.81 it is easy to conclude from the Proposition 13.21 that if 
the data of interpolation problem (j3.18p correspond to a CMV matrix, then either 
Hhrn) = 2m or h{l2m) — 2m — 1. 

Theorem 3.3. Let for the minimal generator r of problem iS.lS^ R^^\0) ^ 
holds. Then l\3.18^ has a unique solution and, hence, the solution of the MISP is 
unique. 

Proof By Theorem [2H] 

A(i)(z) = 5(z)i?(i)(z) + r(z)Q(i)(z), 
A(2)(z) = 5(z)i?(2)(z) + T{z)Q^^\z), 



^Again, we assume that C = C(oo, • ■ • , cin-2', 1) with given oq, . . . , a„—m—2 and the eigenvalues 
Ci > • ■ • : C2m does exist 
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where r and q are the minimal and second generators for the problem (I3.18|) . re- 
spectively. Proposition 13.21 reads that either h{r) = 2m or h{r) = 2m — 1, so by 
Theorem 12. 101 either h{q) ~ 2m + 1 or h{q) — 2m + 2. In the first case 

deg i?(^^=m, degi?(2^ < TO- 1, degQ^^^ < m, deg Q(^^=m, 

and in the second one 

degi?^^^ < m - 1, degi?^^^ TO - 1, deg Q^^^ = m + 1, degQ^^^ < to. 

In view of the degrees of A^^^ and A(^\ in both cases 

deg5(z) = l, degT(z) = 0, 

i.e., for A^-') we have 
(3.22) 



A(i)(z) = {az + b)R'^^\z) + cQ^^Hz), 
A'^^\z) = {az + b)R^^\z) + cQ(^\z), 



For fixed r and q, it is easy to see that a and c are uniquely determined by condition 
(IX^ . If i?(i)(0) ^ 0, then b is also uniquely determined by dMJ), but if Q^^^O) = 
0, then problem has either no solutions, or infinitely many solutions. But since the 
data of the problem is taken from the CMV matrix, the solution does exist, so if 
R^^\0) — 0, the interpolation problem has infinitely many solutions and, hence, 
the MISP may have infinitely many solutions. □ 

In view of this theorem, two natural questions arise. 

1) Is it possible for the minimal generator to have i?'^^^(0) = when (|3.17p is 
related to MISP for a certain CMV matrix? 

2) Is it possible for the MISP to have more then one solution if i?(i)(0) = 0? 

Both answers are positive, so in some special cases MISP with non-unique solu- 
tions does exist, although the number of the "pieces of information" in the inverse 
data is equal to the number of parameters to reconstruct. We provide examples for 
both possible cases /i(/2m) = 2to — 1 (example 1) and /i(/2m) — 2m (example 2). 

Example 1. Let — 1 < 6 < 1 and C = C(0, 0, b; 1). By the Szego recurrences 

$i(z) = z, $2 = z\ $3 = ^3 - b, $*(z) = ~bz^ + 1, 

and 

$4(^) = z^siz) - $*(z) - (z^ - l)(z2 + bz + 1), 
so the eigenvalues are 



Cl.2 = ±l, C3,4 = -2±*yi- J- 

We see that the pair (^i , C2 does not determine b uniquely, although any other pair 
does. 

However, the MISP of (non-unique) reconstruction of b by the two eigenvalues 
Ci^2 = ±1 is still possible. Find the right-hand side of (|3.17p for this case. First, 
consider the Weyl function of the "known" left matrix C(0; 1). Its Szego polynomials 

are $1(2) = z and $2(2) = z$i — 1 • $1 = — 1. So, the reciprocal of its Weyl 

z2 - 1 

function is W{z) — and the right-hand side of (|3.17p is 

z 



. (i + o){i- 0-01^(0)} -(1-0) _ . -QWjQ) 

1-0- omo) ' i-omcr 
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For = ±1 we have ujj — 0, so, according to (|3.17p . the reciprocal of the Weyl 
function of the matrix Cr = C{—b; 1) satisfies 

Wri±l) = 0. 

Remind that it also must satisfy additional conditions ()3.20p and (|3.2f |) with m — 2. 

So, following the procedure of solving the MISP, described above, to reconstruct 
the inverse of the Weyl function of C{—b; 1), we need to reconstruct a rational 

p(i) 

function j^f-^ such that 

p(i)^p(2) are monic and P'-^^O) = 1. 
The corresponding interpolation problem for the vector-functions is 
(3.24) P^^\±1)^0. 

According to Proposition 13. 2[ the minimal generator of the problem p.24p must 
have the height > 1, and according to Theorem l2.81 it must have the height < 2. In 
fact, it can be immediately checked that the non-trivial vector- function of minimal 

height, corresponding to this problem, is K ) , whose height is 1. Further, since the 



^1^ 

height of the minimal generator is 1 and there is 2 points of interpolation, according 
to Theorem l2.10l the second generator must have the height 4 and there is no solu- 
tions of height 2. In fact, it is evident that the vector-polynomial ~^ "'^■jj^ 

whose height is 4, solves (|3.24l) . and there is no solutions of height 2. Finally, the 
general solution of ()3.23p . is 

P^^\z) _ l-{z + l){z-l) + {z + a)-0 _ (z + l)(z-l) 
P(2)(z) ~ i.0-h(z-h6)-l ^ z + b 

with arbitrary number b (cf. p.22p ). However, only those solutions with additional 
condition |&| < 1 give us not only a solution of (|3.23p . but also the Weyl function of 
a CMV matrix of the type C{—b; 1). In fact, let us find directly the Weyl function 
of C(— 6; 1). Its Szego polynomials are: 

Ai(z) ^ z + b; 

A2{z) = zAi(z) - Al{z) = z{z + 6) - {bz + 1) = (z + l)(z - 1). 
So, 

^ Aj ^ + 

Ai z + b 

as was to be checked. 

Example 2. Consider a family of CMV matrices of order 4: C(0, —y, —x; 1); 
— 1 < x,y < 1, and analyze the MISP of reconstruction of the unknown x, y by the 
four eigenvalues. Calculate for them the Szego polynomials: 

$i(z) = z, <fj(z) = l; 

z^i + y^*i = + y, $2 = i + ^^y; 

$3 = z^2+x^2 = z{z'^+y)+x{l+z'^y) = z^+xyz'^+yz+x, $3 = 1+xyz+yz^+xz^; 
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$4(2) = ^ ^ ' ^3 = + xyz'^ + yz + x) — 1 — xyz — yz"^ — xz^ 
= (z^ - 1) + [xy - x)z{z^ - 1) = {z^ - + {xy - x)z + 1). 
Introducing the notation 

(3.25) k :— xy — x, 

we express the eigenvalues of C as 



(3.26) S: Ci,2 = ±l; C3,4 = -^±*yi-^; Ca ^ C4 and C4 = C3. 

Hence, if x and y are related by (13.251) and k is fixed, we have an infinite family 
of CMV matrices C{x,y, 0;1) with the same spectrum S p.26p . According to 
the general theory the auxiliary Weyl functions Wr{z) and W{z) of the matrices 
Cr = C{x,y; 1) and C(l), resp., take the same values on E for different a; and y, 
related by (|3.25p . Find the right-hand side of (|3.17p for this case. First, consider 
the Weyl function of the "known" left matrix C(l). Its Szego polynomials are 
<i>o(2:) = 1 and $i(z) = ~ I ■ — z — 1. So, the inverse of its Weyl function is 
W{z) — z — 1 and the right-hand side of p.l7p is 

= • a + o){i-o-WO)}-(i-o) ^ . -0(0 - 1) .0(1-0). 
' i-o-c,w{z,) 1-0(0-1) 

Let us directly check that the left-hand side of p.l7p coincide with the obtained 
numbers. The Szego polynomials of Cr = C{x, y; 1) are 

Ai(z) = z — X, A^(z) = 1 — xz; 
A2{z) = zAi - yAl = z^ + {xy - x)z - y, ^-^{z) = -yz^ + [xy - x)z + 1; 
A3(z) = ZA2-A2 = z'^+{xy-x)z'^~yz+yz'^~{xy-x)z-l = {z-l){z'^ +{k+y+l)z+l)] 

A3(z) (z-l)(z2 + (fc + y + l)z + l) 



Wr{z) 



K2{z) z^ + kz- 

Wr{l) = = tJi; 



-2(2-(fc + y + l)) 

Wr{-1) = ; = -2 = UJ2] 

l-k-y 

and, since Q + kQ + 1 = 0; j = 3, 4, we have 

M^.(0) = ^^^-§^ = 0(1-0) = -.; .-3,4. 

If we solved the MISP following the procedure described above, we have to find 
(non- uniquely) the inverse Weyl function Wr{z) from its value in the four eigen- 
values and the additional conditions for the numerator and denominator. In this 
example we will restrict ourselves by illustrating that the minimal generator of the 
corresponding interpolation problem does not satisfy the conditions of theorem [ 
which actually cause the existence of infinitely many solutions. 
Consider the interpolation problem, corresponding to this case: 

i?(i)(0)-c^,i?(2H0); J -1,2,3,4; OeS; 

Wl = 0, UJ2^-2, W3,4 = (3,4(1 - (3,4)- 
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(since (j ^ 1, we have ujj ^ 0). We are looking for the solution of height 4: 

rW = + az + (3, + so 

1 + a + /3 = 0; 
1 -a + /3 + 2(-7 + (5) = 0; 
C| + aCj (70 + <5) = 0, J = 3, 4. 

Since + /cC,- + 1 = 0, the last 2 equations can be rewritten as 

(a - k)Q oj,{^C, + <5) = 0, (a - /c)0 - a - 2 - oj,{^Q + <5) = 0; j = 3, 4. 

Exclude from these equations first 7, then 5: 

{a - fc)C3^^4 - (a + 2)w4 - 7C3^^3'^4 - (5w3W4 = 0; ^ 
(a — k)C,j^ijj^ — (a + 2)w3 — 7C4a;4W3 — Suj^ujy^ — 0. 

(a - fc)(C3'^4 - C4'^3) - (a + 2)(tJ4 - W3) - 7|w3p(C3 - C4) = 0. 
(a - fc)(C3 - C4) - (a + 2)(a;4 - ^^3) = 7|1 " ^3p(C3 - C4); 
_ a~k , , 0^3 - u;4 _ a-^ , , , (1 + fc) _ (a + l)(fc + 2) 

^ " |1 - C3P ^ (C3 - C4)|l - C3P " |1 - CsP^^""" V - C3P " |1 - C3P ■ 

But |1 - C3p = 2 - C3 - C3 = 2 + fc ^ 7 = a + 1. 



1) 



2) 



(a - fc)C3 • C4W4 - (a + 2)C4CJ4 - 7C3a;3C4W4 - <5w3 • C4W4 = 0; 
(a - k)Ci ■ C3W3 - (a + 2)C3CJ3 - 7C4W3C4^4 - Suj^ ■ C3W4 = 0. 



(a - fc)(L^4 - C^3) - (a + 2)(C4C^4 - C3C^3) " 5\l - C3|'(C4 - Cs) = 0, 
C4^4-C3C^3 = Cl(l-C4)-C|(l-C3) = (C4-C3)(C4 + C3-C|-C|-1) = (C4-C3) (-^'-^+1), 
C^4 - W3 = (C4 - C3)(l - C3 - C4) = (C4 - C3)(l + k). 

Hence, 

(a - k){\ + k) + (a + 2){k^ + fc - 1) = (5(2 + k); 
ak{k + 2) + (fc - l)(fc + 2) (5(fc + 2) ^ (5 = afc + fc - 1. 
3) We have from the second equation a — S — "f^a = ak + fc — 2 =4> 

/3 = 0; 7 0; S = -1, 
R^^\z) = z^~z- R^^\z) = -l, 
It can be immediately checked that the solution of height 3 (such that i?^^^ — az+(3, 
= z + 7 does not exist, so r = ( j minimal solution. As we see, 

i?(i)(0) =0. 

Finally, rewriting 

(z-l)(z^ + (fc + l)z + l) + (z^-z).y 
•^^^^ z^ + kz+{-l)-y 

we see in the non-determined (arbitrary) up to a constant factor terms of the 
numerator (z^ — z) ■ y and denominator — 1 • y the components of the minimal 
generator (cf. (|3.22p with b = y). 



Remark. Assume that 2m + 1 eigenvalues C17 • ■ • i C2mi Cim+i are known. Since 
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the interpolation problem W2{Cj) — ^j', j = 1,2,..., 2m + 1 has obviously the 
unique solution. 
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